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We carry out the extension of the covariant Ostrogradski method to fermionic field theories. 
Higher-derivative Lagrangians reduce to second order differential ones with one explicit independent 
field for each degree of freedom. 



I. INTRODUCTION 



■ Higher derivative (HD) field theories appear in many physical situations such as the Higgs model regularizations 
[1] or generalized electrodynamics [2]- [4], but their main interest resides in their use as gravitational actions. HD 

. gravitational field models arise as effective low energy theories of the string [5] or are induced by quantum fields 

04 ' in a curved background [6]. Theories with quadratic curvature terms have been studied closely because they are 

(— ( , renormalizable [7] in four dimensions. This property led to Renormalization- Group analyses in [8], which culminate 

^ in [9], including attempts to avoid the appearance of Weyl ghosts usually occurring in HD theories [10]. For these 

, reasons, it is in the gravitational framework where essential progress has taken place. In particular, during the eighties 

■ a mechanism has been devised to deal with the Hamiltonian formulation of an arbitrary HD gauge field theory which 
\ was successfully applied to HD gravity [11]. Somewhat later, in the nineties, a covariant differential order reduction 

OA ■ for HD field theories was found [12]. This kind of covariant techniques were used in the late nineties to identify the 

' propagating degrees of freedom (DOF) in both diff-invariant [13] and gauge fixed HD gravity [14]. Moreover, some 

, effort has been devoted to study general bosonic HD free theories [15, 16] as useful testbeds for HD gravity. 

^Kf^ , HD fermionic field theories have also been considered in the literature , for example in the context of the effective 

• action for the trace anomaly in conformal field theory [17], as a dynamical mechanism for fermionic mass generation, 

OA and in the frame of Faddeev-Popov^ compensating Lagrangian for HD gravity [14]. Therefore it would be very useful 

■ to generalize all the work done in the bosonic case to cover the presence of fermionic fields. This is one of the goals of 
\ the present paper. Another is to provide a general framework for the differential order reduction methods developed 
I for HD theories in [12] and treated in [13]-[15]. The starting point in the study of any field theory, including HD 

O ^ ones, is the characterization of the propagating DOF. There are several ways to do this. The first, and standard, 
D follows a detailed analysis of the Ostrogradski phase space [15], [18] in order to characterize the reduced phase space 
[19], i.e. the subspace of the phase space where the physical degrees of freedom (DOF) reside. However when the 
^ , theory is linear there exists a second and highly useful shortcut provided by the use of the covariant symplectic 
techniques of Witten and Crnkovic [20, 21]. The idea of these methods relies in the construction, directly through 
. , the action, of a symplectic form on the space of solutions to the field equations (covariant phase space) and use it to 
. generate conserved quantities (energy, angular momentum and so on) that characterize the propagating DOF. The 
covariant symplectic techniques have proved to be an essential tool in the classification of free theories and also in the 
identification of the propagating DOF [22]. We will see that they also allow us to complete, with all generality, the 
order reduction in HD theories. 

When we are dealing with HD relativistic field theories, the story does not end once the DOF have been identified. 
The reason is that HD theories can be usually reinterpreted in terms of lower derivative (LD) ones that propagate (in 
the free limit) according to standard Lagrangians. The machinery necessary for this reinterpretation combines the 
use of the covariant Legendre transform [12] and a subsequent diagonalization in the fields [13]. This procedure has 
been developed in various examples: diff-invariant HD gravity [13], gauge fixed HD electromagnetism [4], gauge-fixed 
gravity [14], and HD scalars [15]. However, there are at least two issues left aside in these works. The first one is the 
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^ Fermionic theories in which the field equations are differential equations of order greater than one are referred as HD fremionic theories. 
Although the Faddeev-Popov compensating fields for HD-gravity are not femionic fields (their spins are or 1) they are anticommuting 
ones. In this sense they are different from the usual bosonic fields treated in [15, 16]. 
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extension of the above results to ferniions, covering the general derivative case^. The second and more interesting is 
to find a unifying point of view that encompasses all the previous works. The motivation of the present paper is to 
clarify these two points. 

The paper is organized as follows. After this introduction, section II is devoted to review the use of the covariant 
symplectic techniques for Dirac fields and to fix some conventions. Section III deals with a simple HD fermionic field 
theory in order to make it simpler to understand the general iV-derivative fermionic theory -under the hypothesis of 
non-degenerate masses- treated in section IV. We end the paper with several comments and our conclusions (section 
V). Some details of the computations and the general conventions are left to the appendices A and B. Finally 
appendix C indicates how the differential order reduction method proposed in section IV can be extended to cover a 
wide class of HD field theories. 



II. A REVIEW OF THE DIRAC LAGRANGIAN 



Let us start with a quick review of the standard Dirac Lagrangian in order to fix conventions and notations. As is 



well known, the free propagation of spin 1/2 and mass m modes can be described through the Dirac Lagrangian C 
that can be written in two equivalent forms 



D 

ml 
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V'(« p — m)ip — V'(j P + 'm)'^} = ipi'i P — fn)^} — 



where m is a mass parameter, '0 is a Dirac spinor, p := 7'*9^, 7'^ are the Dirac matrices, ip := V^7°) and \ denotes 
complex conjugation and transposition (sec appendix A for a resume of the conventions). The first and more symmetric 
expression is appropriate for analytical purposes and will be used in the following sections when we introduce the 
differential reduction order methods. The second is suitable to define the propagator. The Euler-Lagrange equations 
associated to the Dirac Lagrangian are 

{iP-m)tp = Q, (1) 

and the oriented ■(/'-^ fermionic propagator is 

1 id -\-m 



m □ + w? 



The dynamics is governed by the field equations (1) that define the Dirac covariant phase space S^. To parametrize 
<S^ we solve the field equations using spatial Fourier transform and the notation introduced in appendix A. It is 
straightforward [23] to prove that ip G 5° if and only if 

V'C*.^) = E Xa K(fc)n«(fc)e-^'=- +6«(fc)z;«(fc)e''=-] , (2) 

where u°'{k), v"{k) are basic spinors, kx := k^t — k-x, :— \/k'^ + vn? , and the functional parameters that distinguish 
one solution from another are encoded on the a"(fc) and b°'{k) fields. 

Let us show how the symplectic covariant techniques [20] induces a symplectic form fl^ in the space S^. First 

define the 2-form 



'E 

where da^ is the measure on a space- like hypersurface S, 



dC" — dC" — 



^ Within the frame of HD scalar theories and covariant Legendre transform, reference [15] fails to cover the diagonalization of the general 
derivative situation. The problem has been solved in [16] avoiding the covariant Legendre transform by using Lagrange multipliers and 
symplectic covariant techniques. 
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and (f , A are, respectively, the exterior derivative and the wedge product in the linear manifold 5"^^. 

The density is divergence- free when restricted to iS^, in other words d^cuo = modulo field equations. Then 
the restriction 0° of 0° on is a well defined 2-form on this functional space, i.e. is time-independent, and in the 
parametrization (2) can be written as^ 

^™ = ^ E / ^^Ua"*ik)Ada"ik) + clb^k)Adb-*ik)]. (3) 

The symplectic form (3) determines the propagating DOF of the Dirac theory -the canonical pairs a"*(fc)-a"(A;) 

and b"*(k)-h"{k) and can also be used to derive constants of motion if one considers the following version of the 
Noether theorem. When the symplectic form f2 is invariant under a group of transformations and we take a vector 
V tangent to an orbit of this group it is straightforward to prove [21], that locally zyil = dH, where iyil denotes 
the contraction of V and il. The quantity H is the generator of the symmetry transformation corresponding to 
V. If the action is Poincare invariant we obtain in this way the energy-momentum and the angular momentum 
densities (with the right symmetries in their tensor indices) by computing iyfi for vectors V describing translations 
and Lorentz transformations and writing the result as dH. Noticing that under a translation of parameter we 
have ivTda°'{k) = iT^kij,a°'{k) and ivTdb°'{k) = iT'^kfj,b°'{k), the energy of a solution is given by 



/ cf fC 7T7 r — ♦ — ♦ — * — ♦ 



It is worthwhile to point out that the minus sign between the a and the b fields forces us to take them as anti- commuting 
variables to ensure de positivity of the energy. 

III. A SIMPLE HD FERMIONIC THEORY 

Once we feel comfortable with the notation, lets us consider the following simple HD fermionic Lagrangian 

^m\m2 = V'(*T - rni){if - m2)tp , 
with real mass parameters^ mi < m2- It is straightforward to write down the HD-field equations for this model 

{ip - nii){iP - m2)ip = , (4) 
whose space of solutions, S^^^^, can be parametrized in terms of a sum of Dirac fields in the form 

= Vi + V'2 , (5) 

where {ip — mi)tpi = 0, for I = 1,2. Explicitly, in the notation given in Section II, 



where, as in previous section, vf^k) are basic spinors defined in appendix A, kix = kit — k-x, kf := yk"^ + mf , 

and the parameters that distinguish one solution from another are encoded on the af{k) and bf{k) fields. Also the 
HD propagator can be easily found as a sum of Dirac propagators 

A'"' = = f A° - A° ) . 

'"1"*= {ip - mi){ip - m2) m2-mi^ 

The physical interpretation is clear. The theory describes two LD fermionic DOF, a physical one (positive contribution 
to the energy) with mass m2 and a Weyl ghost (negative contribution to the energy) with mass mi. 



* z* denotes the complex conjugate of z. 

^ The problems associated with the presence of degenerate and/or complex "masses" are detailed in reference [16]. 
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We can learn more about this theory by looking it through the glass provided by the covariant symplectic techniques. 
The extension of the covariant symplectic techniques to the HD field theories was done in [24] . Following [24] the HD 
Lagrangian has an associated 2-form 



where 



= -i{mi + m2)(Ii^ Ajf" dTp + {d'' dip) Ad^p - dip Ad^ dip. 

It is easy to show that lu'^^) is real and, modulo field equations, S^wJ^j = 0. Then we can compute the symplectic form, 
^mirn2' Space of propagating DOF, S^^\^^, making use of the parametrization given by equations (5)- (6). 

This leads to a well defined 2-form on S!^^^^, namely 

and Noether's theorem gives us the energy 

Hi^l^^ = {m2 - mi) {H^^ - H^J . 

This supports the intuitive interpretation of ^mim2 S^^^t^ by the propagator decomposition as a theory describing 
two Dirac fields. In fact the simplectic form confirms that this is the right interpretation and also, through the 
Noether theorem, that the field ipi is a Weyl ghost. It contributes to the energy with a wrong sign, destroying the 
semi-boundedness of the Hamiltonian and consequently the unitarity of the quantum formulation of the theory. 

Finally let us show how it is possible to find a covariant Legendre transform connecting the HD theory with a 
LD one where the usual Dirac fields are explicit. To this end we define the first order differential operator D := ip, 
that satisfies ipi{Dip2) = J^4 {Dipi)ip2, and rewrite the Lagrangian in terms of this object. Namely, modulo total 
derivatives, 

(V'. Di,, Dip) := DlpD^ - Vll^^ + + mim2?^V • 

Then we introduce a generalized Legendre transform [12] with respect to the momenta generated by the D-operator^ 

a^v 2 ^ ' dBti) ^ 2 ^ ^ ' 

This transformation is not singular and it permits the inversion of the derivative of as a function £>'^ = z/(7r, ipi) 
given by 

mi + m2 — -T- _ mi + TO2-r 

j/(7r, V) := TT H ^ Ip ; v{'k, ■(/)):= tt H ^ tp . 

The Legendre transform (7) has an associated Hamiltonian 



'^mim2(V'.V',7r,7r) := Trv{'K,ip) + iy{n,ip)'!r - i:^^^„^{tp,ip,u{TT,2p),v{Tr,ip)) 

_ mi+m2 ,-r , _ , X / mi - ma y - 
= ttttH [ipn + 'Kip) + [ I tpip , 

and the HD dynamics (4) can be described now by means of the Euler-Lagrange equations derived from the Helmholtz 
LD Lagrangian 

i-iP, Ip, TT, 7f ; D^p, Dip) := WD^p + Dlpn - H'Zm, (V', ^r, 7f) . 



We must treat a field and its conjugate as independent variables in the Legendre transform in order to allow the presence of grassmanian 
fields [14]. Then, the derivatives with respect D'tji and Dip are respectively left and right derivatives. 
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Although L'^^ is a LD Lagrangian classically equivalent to i^J^'^^^i it does not explicitly exhibit the propagating DOF. 
However the diagonalization can be carried out taking into account the special structure of the covariant phase space. 
Defining the new fields '^i and V'2 in the form 

i) =: ^i+V2, (8) 

TOl - 7712 , , , , 

TT =: 2 (^1+^2) , (9) 

the Helmholtz Lagrangian decouples the propagating modes 

£g' = (m2-mi)(C^-Cj . 

The relative minus sign shows the presence of a Weyl ghost. 

The diagonalization (8)- (9) can be found by means of a very natural reasoning. Over the solution space, where the 
DOF reside, we can write V = "01 + V'2 with (i^ — m;)^'; = 0. Then, on this space 

/ / , / n / rni+m2 , rm - m2 

V' = Vi+V'2 ; 7r = i»V ^ — ■0 = — ^ — (V'i-V'2)- 

The above relations, that coincide with the proposed diagonalization, indicates how the DOF are encoded within the 
fields tjj and tt. Moreover, the relation {■^,7r} <-> {V'ij^i'2} is invertible because of the invertibility of the Legredre 
transform (7). 

We are ready to generalize the above results, obtained in a very simple 2-derivative framework, to the general 
iV-derivative situation. 



IV. iV-DERIVATIVE FERMIONIC THEORY 



In this section we consider the general A''-derivative Lagrangian 

N 

1=1 

where {ip — mi) are Dirac operators with real mass parameters m; that have been ordered following the Z-index 
ordering, i.e. m; < m;/ when I <l' . 

The propagating DOF described by can be read out from the algebraic decomposition for the HD propagator 
in terms of Dirac propagators 

1 ^ A" 

Notice that the sign alternates in the coefficients Wi^iirm — mi') so the occurrence of Weyl ghosts is to be expected. 
Mathematically the propagating DOF are points in the space of solutions to the field equations that we refer in the 
following as The covariant phase space is defined by the Euler-Lagrange equations 

N 

]J(z^-m,)V = 0. 
1=1 

These equations can be solved by means of a linear combination of Dirac fields in the form 

N 
1=1 

where 

(ip - mi)iJi = . (10) 
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As in previous sections, the standard parametrization for the ipi fields satisfying equation (10), that provides us with 
a exphcit parametrization of , is 

Mt,x)=J2 I 7|^5[ar(fcKWe-*^ + 6r(fcK(fc)e*-] . (11) 



The decomposition of the HD propagator as a sum of Dirac propagators, with alternating sign coefficients, and 
the decomposition of the space of solutions as a direct sum of Dirac spaces indicates that the theory represents the 

propagation of A^-Dirac fields, some of them physical and some of them Wcyl ghosts. This is the case, and it is possible 
to give a more precise proof of this fact. The Lagrangian induces a symplectic form in the space of spinor fields 
by means of 



where 



/firm r)r<~> r)r("> 



+ •••+(/ jf, A(/i/v,....^_, 

This is so because over the <S^> space the 2-form density wj^) satisfies d^uj^^^ = 0. Hence, the restriction of O'"^ to 
'^m ' that we refer as fi*^' is a well defined 2-form on this functional space, i.e. is time-independent. In fact, it is 
straightforward, but highly tedious (see appendix B), to compute this restriction in the parametrization given by (10) 
to obtain 

1=1 Ui^i'i'nT-i-'n^i') 

Consequently the DOF of the theory are a sum of Dirac DOF. Finally the energy, computed through Noether theorem, 
takes the form 



Hi::' = E 



N TTD 

mi 



Because of the sign alternates in the coefficients ri(/('(™' ~ namely sg ^ri^(/(TOi — tny^ = (—1)^+', some of 

the fields are physical (positive contribution to the energy) and some are Weyl ghosts (negative contribution to the 
energy). 



Once we have identified the propagating DOF -and due to the special form of the propagator, the covariant phase 
space, the symplectic form, and energy that can be obtained as linear combinations of Dirac objects- it is plausible 
to presume the existence of a mapping that transforms the original HD theory into a sum of LD Dirac theories. To 
this end, as we did in section HI, it is convenient to follow a series of preliminary steps. First define the differential 
operator D := and then expand the differential kernel n;Li(-^ ~ '^O appearing in in the form 

N N 

l[{D-mi)^J2^iD^~'' 

1=1 1=0 

where 

Co := 1 ; q := (-1)' ^ ma,---ma,, l = l,...,N. 

ai<---<ai 
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A last technical remark is still in order. The heavy algebra involved in the following forces us to treat the N odd and 
N even cases separately. We work out in detail the N = 2n one and refer to the bosonic framework considered in [15] 
and [25] to understand certain peculiarities present for N = 2n— iJ 

Modulo total derivatives the HD Lagrangian can be rewritten in the more convenient form 



/=0 1=1 

This expression suggests the introduction of the Ostrogradski-like variables 



Xr:=D''-'i, ■ x^:=£»-i^ , r = l,...,n, (14) 
and their corresponding momenta by means of a covariant Legendre transform 

ci _ dC Cl_ 

DXn + -TTXn ; 7r„ := — — = DXn + irX„ • (15) 



„ _ dC ^ , , 

It is easy to show that this Legendre transform is non-singular due to the invertibility of the highest derivatives Dxn 
and Dxn in terms of the x and tt variables, namely 

l^niX,'^) ■= DXn = T^n-^Xn l'n{X,'^) ^Xn = - ^X^ ■ 

Thus it is possible to find an expression for the Lagrangian in terms of this new variables 



n-l n-2 

'Cm"^ = T^nT^n - -^XnXn + X] C2TXn-rXn-r + ^ '^n-rXn-l-r + Xn-l-rXn-r\ 

r=0 r=0 



The Ostrogradski-like Hamiltonian associated with the Legendre transform is 

n — 1 n — 1 



"Wm"' := 7r„f„(x, tt) + Vn{X, '^)'^n + "^rXr+l + ^ Xr+lT^r - jO.^^\X, X, ^niX, T^) , ^n{X, T^)) 



r—1 r=l 

— Cl _ n C?_ 

= 7r„7r„ - — [KnXn + X„7r„J + —XnXn + 

n — 1 n— 1 n— 2 

^"irXn—rXn—r / ^ „ [Xn—rXn—l—r + Xn—l—rXn—r\ ■ 
r=l r=0 r=0 

Finally the Helmholtz Lagrangian -the LD Lagrangian, clasically equivalent to C^^\ generated by the Legendre 
transform- is 



(2„) 

m 



As usual in the covariant Legrendre procedure [15, 16], the expression for does not exhibit the propagating DOF. 
However, its straightforward to find a new set of variables in terms of which the propagation is explicit. Specifically 



In the AT-odd case the definition of the highest momentum yields a constraint -that can be easily taking into account- while the field 
derivative is worked out from the next momentum definition. Once the constraint is solved, the odd case is exactly analogue to the even 
one. 
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it suffices to introduce the new set of ipi-&elds^ through 

2n 

E^rV/, (16) 



Xr 



1=1 

2n 

Y^^r'i^i, (17) 

2n /2(n-r) 

2n-r-k , '^2(n-r) + l r-1 



TTj. 



= E E -fc-?"-^-' + -^^^mr' , (18) 

1=1 y k=o J 

2n ( 2{n-r) \ 

= E E + "-^^^rnr i^i . (19) 



2n-r-k I ^2(n-r) + l r-1 \ "T 

cfe/; 

1=1 \ k=0 

In this variables, and modulo total derivatives, the Helmholtz Lagrangian becomes diagonal 

2n 2n 

1=1 1'jil 



The key observation to find the diagonalization (16)-(19) relies in the following simple fact. The Ostrogradski variables 
Xr and momenta tt,. can be expressed by means of a one-to-one map in terms of the original HD variable ijj and its 
derivatives D'^ip, namely 

Xr = D^-'i;, (20) 

2(n-r) 

TT, = E CfcZ)2n-r-fc^ + £H(!iZ!:)±ii)r-l^_ (21) 
fe=0 

When we restrict the relations (20)- (21) to the space of solutions parametrized by (11) -that is, over the propagating 
DOF space- we reobtain the one-to-one linear relations (16)-(19) between the {xr^T^r} and {-0;} variables. This is 
the reason why this linear redefinition diagonalizes Cf^. 

All the preceding results can be generalized to cover a wide class of linear theories including those considered 
in previous works [4] and [16]. In particular the procedure followed here resolves the deficiencies inherent to more 
primitive approaches [15]. Specifically it leads us to the diagonalization for the Af-derivative theory within the 
framework of the covariant Legendre mapping. We summarize in appendix C the essential steps and requirements 
that allow us, following the lines of this section, to transform a HD linear theory into a LD one where the propagating 
DOF are explicit. 



V. CONCLUSIONS AND COMMENTS 



We have proved the equivalence of HD fermionic field theories and a LD counterpart where the DOF are explicit 
by means of a covariant Legendre transform [12] and a subsequent diagonalization. The previous attempts to solve 
this kind of problems, that considered only scalar field theories [15], failed to find the general diagonalization to an 
arbitrary differential order due to the heavy algebra involved. A way out of this was given in [16], in a more general 
framework than [15] -and also dealing with bosonic field theories-, but the solution proposed there abandons the use 
of the Legendre transform in favor of the Lagrange multipliers method. The use of the covariant symplectic techniques 
in combination with the covariant Legendre transform avoids the algebraic problems of [15] permitting us to find an 
explicit formula for the DOF diagonalization at every differential order and to generalize the previous results to the 
fermionic case. 



At this point the fields are simply a new set of variables. We follow the same notation used to denote the Dirac fields in (11). We 
will show that these new fields diagonalize the Helmholtz Lagrangian in terms of the Dirac ones. 



9 



The approach that we follow here is exportable to more general theories -classical mechanics, HD gravity, and so 
on- as we schematically show in appendix C. It is also important to realize that the presence of gauge symmetries 
in the theory docs not change the results in any way. When dealing with HD theories written in terms of differential 
forms, such as those considered in [16] , the role of the I^-operator {D = ip in the present work) can be simulated by 
an operator constructed by means of the exterior differential d and its dual 5, for example D = 5d. Even considering 
the possible existence of gauge symmetries {DdA = 0, for any field A) the analogue of the Legendre mapping (14)-(15) 
is still nonsingular and therefore the reduction order procedure works exactly as in section IV. 

Another remarkable fact is that, in spite of the fact that it is always possible to reduce the differential order by 
means of a covariant Legendre mapping, the diagonalization of the resulting LD theory in terms of a sum of standard 
theories it is not always possible [16] (but when possible, the method presented here allow us to find it). To fix ideas, 
let us return to the fermionic case where the more general HD, Lorentz invariant, Lagrangian for a spinor tp without 
internal indices takes the form 

=^{D'^+ CiZ?"-l + • • • + Cn-lD + Cn) i> , 

where Ci are real parameters with the appropriate dimensions. As is well known, the polynomial + cil?"^^ + • • ■ + 
Cn-iD + Cn can be factorized in terms of its roots. In section IV we have make the assumption that all the roots are 

different^. This is a necessary hypothesis in order to succeed in the diagonalization process. However, if we consider 
a HD "mass-degenerate" model, the diagonalization cannot be carried out. For example, starting with 

the Legendre procedure presented in section IV leads us to the LD Hemholtz Lagrangian 

^s-(?w)(,!„--'»)(t). 

The propagator associated with C^^ is 



I D-m 



{D-my\D-m 



As it is clear, this propagator cannot be diagonalized by means of a linear redefinition on the fields (V', tt) not involving 
differential operators. This is so because l/(-D — m)^ is not a linear combination of Dirac propagators and then it 
is not possible to find a c-number matrix Q such that (5^A(^7r)Q becomes diagonal with Dirac propagators on its 
diagonal elements. The same conclusion can be reached with an elementary analysis of the covariant phase space. 

Finally, we want to point out that the diagonalization process of the LD equivalent theory, obtained by Legendre 
transform, strongly relies on the decomposition of the HD solution space as a direct sum of LD ones. This decompo- 
sition is tied to the decomposition of the HD propagator as a sum of simpler pieces or, at least, on its relationships 
with simpler propagators. Then, in order to make practical use of the differential order reductions presented in the 
paper, the knowledge of the properties of the Green functions of HD differential operators is needed. Some results in 
this direction can be found in [26] where, by means of heat kernel methods, the Green functions of a wide class of 
products of second order differential operators have been studied and their relationship with LD counterparts pointed 
out. These results are a key ingredient if one wants to find the LD-diagonalized equivalent of such a HD theory by 
means of a covariant Legendre transform. 

APPENDIX A: CONVENTIONS 

We use the Minkowski metric rj^u = diag{+l, — 1,— 1,— 1) to lower and raise space-time indices and use the Ein- 
stein convention for summation over repeated indices. The derivatives with respect to coordinates are sometimes 
abbreviated as 



® Wc have assumed also their positivity. This is a necessary requirement if we want these parameters to be physical masses, but it can be 
relaxed for the diagonalization purposes. 
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and □ := d^d^^ =8^—8^ is the d'alambertian operator. 

For any X'^'^'"^'', X^'^^'"'^''^ represents its index symmetrization, namely 



k\ 

Trent'') 

where n'*"' denotes the permutation group of fc-elements. 
The Dirac 7 matrices satisfy 

with 7° hermitian and 7' antihermitian. As usual ^ := k^'y'^. 

Dirac spinors u°'{k) and v°'{k) are a basis of solutions to the Dirac equations 

{^-m)u''{k) = ; {ft + m)v"{k) = . 

They depend implicitly on the mass m and explicitly on the on-shell momentum k with + m^, and the a 

index that labels the polarizations. In this basis, and making use of the spatial Fourier transform, any spinor field 
can be expressed as 

d^k m 



i,{t,x)=Y. I 47\^'^W{k,t)u'^{ky^-^ + h'^{k,t)v'^{k)e- 



a=l,2- 

The functions a°'{k,t) and take the special form a"(/c,i) = a"(A;)e-''=°* and = 6"(A;)e''=°* when 

considering spinor fields that satisfy the Dirac equation {ip — m)'ip = 0. 

The basic spinors u and v satisfies the following normalization properties 

u"(A;)u^(fc) = Socp ; v''{k)v^{k) = -6^0 ; v"{k)u'^{k) = u°'{k)v'^{k) = . 

m m 

v"\k°,-k)u^\k\k)^u"\k^,-k)vf\k^,k)^0 , 

where the conjugate spinors u are defined as u := u^7°, and f denotes transposition and complex conjugation. 

APPENDIX B: SOME REMARKS ON THE COMPUTATION OF THE SYMPLECTIC FORM 

Wc summarize here the main steps involved in the calculation of the symplectic form (13) of Section IV. The 
symplectic form is given [24] through space integration of the density 

■■= i{^-d,,^+---+{-ir-'8,,...d,^_,-^ — A</^ 

^^W^- + • • • + • • • 8ir^ — ) ^'^^^ 

X'-'YUvi 'JynviiH ^YiJ,i>ini---iJ,N-2 / 

+ •••+ d— ^ AcfV'.,....„_i . 



In order to compute the derivatives of the Lagrangian that appear in the definition of wf^j is convenient to rewrite 

N N 

^P^,,..., := i^U^iP - m)i> = '^N-k^W^ , 



fe=o 
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where the constants Cfc are defined in terms of mass products 

Co := 1 ; Cfe := (-1)*= ^ ma^---ma 

ai<---<ak 

Thus, in this notation, 



c;v-fe^W^•••7'"=^ (Bl) 

where the symmetrization of the Dirac gamma products can be expressed in terms of the inverse of the Minkowskian 
metric ry'*'' in one of the following forms 

^{w . . . ji^2k} = ^{piM2 . . . ^iJ.2k-m2k} (even) , (B2) 

/y{^*l . . .,yA'2fc-l} _ ^{ftl^ft2H3 . . . jj^*2fe-2 A»2fc- 1 } (odd) , k GN. 

Phigging (B2) into (Bl) and taking care of the combinatorics, we found that the terms of the density Wf^j belong to 
one of the following four categories 

rir(^) _ 
1. d,,...,,,_, -—^ = (-l)^C;v-2fetf-^a''V . 

0''i/^jUjUi---jU2Jc-l 

Now making use of the parametrization (11), that is V = J2fLo where 

d^k 



tl>i{t,x) 



E / 7^S[<W<We-*" + &rW(fc)e*1 , k? = \/k^ + m^ 



remembering that the constants Cfe are defined in terms of the masses mi and using the properties of the spinors ui 
and vi detailed in Appendix A it is straightforward, but slightly tedious, to derive equation (13) for the symplectic 
form. 

APPENDIX C: THE UNDERLYING IDEA 

The procedure developed in main body of the paper for spinorial fields can be straightforwardly extended to cover 
all linear models for which the field equations can be derived from a variational principle of the form 

S^M = {^\{D + Mi)---{D + M^)4>) , (CI) 

where Mj are "mass" parameters (usually Mj = — rrij when dealing with fermions, Mj = mf when dealing with 
bosons), ordered by their subscript M\ < M2 < ■ ■ ■ < Mn, {■]■) is a pseudo-scalar product, and D a differential 
operator symmetric under (■!•), that is 

(.^1 I Dcj>2) = {D(j>t I (j>2) . 

Gcncrically, in the fcrmionic case D will be a first order differential operator such as D = and in the bosonic case 
a second order one such as D — Sd for differential forms or D = □ (^P*-^-* — P^"') for HD-gravity [16]. 

Under this hypothesis, and through the lines presented in Section IV, it is a trivial task to define a correspondence 
between Sho[(P] and the action 

TV N 

J2Y[iMb-Ma)SM, (C2) 

0=1 bjta 



12 



where 

Sa[(l>a] ■■= {cl>a\{D + M„)(j>a). 

This can be done by following the same steps as in the spinorial case. First define a Legendre transformation of the 

form (14)-(15) and then a linear redefinition in the form of (16)-(19) to select the propagating DOF. Formally, the 
same formulas are valid in the general case with the obvious identifications ip (j), ip and m; —Mi. 
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